Measurement of the magnetic penetration depth of a superconducting MgB 2 thin film 

with a large intraband diffusivity 
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We report the temperature dependent magnetic penetration depth A(T) and the superconducting 
critical field H C 2(T) in a 500- nm MgB2 film. Our analysis of the experimental results takes into 
account the two gap nature of the superconducting state and indicates larger intraband diffusiv- 
ity in the three-dimensional (3D) n band compared to that in the two-dimensional (2D) a band. 
Direct comparison of our results with those reported previously for single crystals indicates that 
larger intraband scattering in the 3D tt band leads to an increase of A. We calculated A and the 
thermodynamic critical field H c ~2000 Oe employing the gap equations for two-band superconduc- 
tors. Good agreement between the measured and calculated A value indicates the two independent 
measurements, such as magnetic force microscopy and transport, provide a venue for investigating 
superconducting properties in multi-band superconductors. 



I. INTRODUCTION 

During the last decade a significant effort has been 
made to understand the mechanism of two-band super- 
conductivity in MgB2ii~— MgB2 has two s-wave gaps re- 
siding on four different disconnected Fermi surface (FS) 
sheets: two axial quasi two-dimensional (2D) u-band 
sheets and two contorted three-dimensional (3D) 7r-band 
sheets. The a band forms two concentric cylindrical 
sheets via in-plane sp 2 hybridization of the boron valence 
electrons. The tt band results from the strongly cou- 
pled covalent bonding and antibonding of the boron P z 
orbitals<£ Multiple bands allow for both inter- and intra- 
band scattering. It is thus possible to tune the upper 
critical field (H c2 ) via doping, which has different effects 
on the inter- and intra-band scattering strengths 
In MgB2 the anisotropy of the temperature dependent 
penetration depth A, j\(T) = A C (T)/A &(T) shows re- 
markably different behavior compared to that of H C 2, 
j Hc2 (T) = H$(T)/HZ 2 (T)M2. This difference indi- 
cates that the two-band nature of superconductivity 
profoundly alters the superconducting properties com- 
pared to those in a single band material^ For exam- 
ple, the equations for critical fields and depairing cur- 
rent as a function of A and £ should be modified due 
to the inter/intra-band scattering. Knowledge of the 
absolute values of A and £ is also important for tech- 
nological applications^ For example, the acceleration 
field in superconducting radio frequency (SRF) cavities 
could be enhanced by covering conventional supercon- 
ducting Nb cavities with superconductor/insulator mul- 
tilayers (such as MgB2) with higher thermodynamic crit- 
ical field (H C )H 

A number of measurements have been performed to 
determine the absolute value of A in MgB2.— The re- 
ported values of A range from 40 nm to 200 nm, indicat- 
ing that A is strongly affected by inter- and intra-band 



scattering In this paper we present measurements of 
the absolute values of A(T), employing low temperature 
magnetic force microscopy (MFM), and of the angular- 
dependent H C 2(T,9) performed via electrical transport, 
in a 500-nm thick MgB2 film. Our MgB2 film can be de- 
scribed by the dirty limit two-band Usadcl equations. We 
analyze the measured values of H C 2 and A using a model 
developed for dirty superconductors^ which simplifies 
the analysis compared to that reported in Ref. 15. We 
investigate theoretically the influence of the inter/intra- 
band scattering on the superconducting properties. Us- 
ing a two-band superconductor model with parameters 
obtained from a fit to H C 2(T,9), we calculate A and H c 
which are consistent with the experimental values. 



II. EXPERIMENT 

A MgB2 film was grown on a r-sapphire substrate by 
a reactive evaporation technique The film is epitax- 
ial and shows columnar growth morphology, with the c 
axis tilted by a few degrees from the normal direction 
of the substrate. For more details see Ref. 22. The 
sample has dimensions L=4 mm x W=5 mm x t=500 
nm, and exhibits a full superconducting volume fraction 
based on measurements using a commercial SQUID mag- 
netometer (Quantum Design magnetic property measure- 
ment system, MPMS) All MFM measurements described 
here were performed in a home-built low temperature 
MFM apparatus^ 3 - Temperature dependent vortex im- 
ages were taken in the frequency-modulated mode after 
a small magnetic field was applied above T c (field-cooled). 
We used high resolution SSS-QMFMR cantilevers^ The 
magnetic field was always applied perpendicular to the 
film surface and parallel to the MFM tip. The absolute 
values of A(T) were determined by comparing the Meiss- 
ner response curves with those for a reference sample at 
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4 Ki 25 i 26 The Meissner technique for the A measurement 
was first proposed by Xu et alwL and demonstrated by Lu 
et al— The film thickness of 500 nm is larger than A « 
200 nm, which makes corrections to A due to the sample 
thickness insignificant. Conventional four- lead resistivity 
measurements used for determining H C 2(T, 9), where 9 is 
an angle between the applied magnetic H and the crystal- 
lographic c axis, were performed with a rotatable probe 
in a commercial Quantum Design physical property mea- 
surement system (PPMS), in magnetic fields between 
T and 9 T. The superconducting critical temperature T c 
= 38.3 K (zero resistance) and the transition width AT C 
= 0.5 K were determined from the transport measure- 
ments. Zero-field-cooling measurements at the MPMS 
with H sal Oe show T c =38.0 K. The small value of resid- 
ual resistivity ratio (RRRw4) indicates the presence of 
impurities, consistent with the dirty limit. 



III. RESULTS AND DISCUSSION 
A. MFM measurements in the MgB2 film 

Figure QJa) presents a typical vortex image in the 
MgB2 thin film. The well-formed vortices in the 6 /nm x 
6 /im field of view were observed, which suggests the ho- 
mogeneity of the sample on a micron scale. However, 
the irregular shape of individual vortices suggests the 
presence of inhomogeneity in the superfluid density on 
a sub-micron scale, which may be related to impurities. 
Figures [U^b) and (c) show MFM images of isolated vor- 
tices in MgPj2 at 4 K and 15 K, respectively. The fea- 
tures besides a single vortex represent a sub- micron scale 
inhomogeneity, indicating small variations of superfluid 
density. Figure HJd) depicts a line profile taken along the 
dotted line in Figs. [TJb) and (c) for each of the vortices. 
The maximum force gradient [max(df /dz)\ at the center 
of the vortex qualitatively indicates that the magnitude 
of A at 15 K is larger than that at 4 Kj22t— In order to 
determine the absolute value of A, we performed Meiss- 
ner experiments. The force between the tip magnetic 
moment (a distance d above the sample) and the shield- 
ing currents induced by the tip field is equal to the force 
between the real tip and the image tip, with the mirror 
plane at a distance A below the sample's surfaced This 
force therefore is a function of d + A when d 3> A. Direct 
comparison of the Meissner curves taken at 4 K for MgPj2 
and a reference sample (Nb) with a known A gives A(4 K) 
= 200 ± 30 nm for MgB2<2£ Comparing Meissner curves 
for MgB2 at 4 K and at a given temperature T yields 
<5A(T). We obtain the absolute value of the temperature 
dependent A(T) by adding SX(T) to A(4 K). Figure 0a) 
shows the Meissner force response as a function of the 
tip-sample distance at several temperatures. The sys- 
tematic evolution of the Meissner response with respect 
to temperature reflects the change of A with temperature. 
Figure &b) shows the normalized A(T) (black squares) 
obtained for MgB2 using the procedure outlined above, 
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FIG. 1: (Color online) (a) A typical vortex image with a tip- 
lift height of 300 nm in the MgB 2 thin film, (b) and (c) Single 
vortex images with a tip-lift height of 300 nm, acquired at T 
= 4 K and T = 15 K, respectively, (d) The single vortex 
profile along the dotted lines in (b) and (c). Higher peak 
value corresponds to a smaller A value. 



deviating significantly from the BCS theory curve (the 
red dashed line), which is consistent with the previous 
studies shown as green solid circles*^ This discrepancy 
indicates a profound effect of two-band superconductivity 
in MgB2 iA£ The large A in MgB2 may be due to inclusion 
of impurities, such as C, N, and Al, which significantly 
affects the electron mean- free path in each band of MgB2 . 



B. H C 2 measurements in the MgB2 film 

In order to investigate the nature of disorder, we per- 
formed temperature dependent H C 2 measurements. Fig- 
ure [3][a) shows H c2 (T) with the field parallel to the c 
axis W^(T) (black circles). The value of H C 2(0) is con- 



siderably higher than that found in clean single crystals 
{H^(0) « 3-5 T )2&, which indicates that the film is 
in the dirty limit. The Gurevich model for two-band 
superconductors 2 ^ considers inter- and intra-band scatter- 
ing by non-magnetic impurities in the dirty limit. The 
high T c in our film (which shows essentially no suppres- 
sion compared to the clean crystals) is consistent with a 
small inter-band scattering, so we can use the equations 
obtained for 7J c2 (T) neglecting the inter-band scattering: 



a 2 [ln(i) + U(r)h)] + ai[ln(t) + U(h)) 
haopn(f) + U(h)][\n(t) + U{r)h)} = 0, 



(1) 



where U(x) = V(V 2 
gamma function, ai = 



- V(V 2 ); ^( x ) is the di " 
L-/L , a 2 = 1 - L-/L , 
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FIG. 2: (Color online) (a) Temperature dependence of the 
Meissner response in MgB2. (b) A(T) marked by the black 
squares are inferred from the data shown in (a). The blue solid 
curve shows the calculated A(T) from the gap equations for 
two-band superconductors. The red dashed curve represents 
the conventional BCS model. The green circles are taken from 
tunnel diode resonator measurements (Ref. 10). 



2w/L , L = J (L 2 _ + 4L 12 L 21 ), L± = L n ± L 22 , 

w = LnL 2 2 - L 12 L 21 , t = T/T c , rj = D 2 /Di, and 
h = H c2 Di/2§ T. $o is a single magnetic flux quan- 
tum, and D\ and D 2 are the intraband diffusivities. The 
angular-dependent diffusivities Di(9) and D 2 (9) for both 
bands are calculated using the following equation: 



D m {6) 



£>L o)2 cos 2 



(2) 



The many body effects such as mass renormalization 
and impurity scattering are encoded in the diffusion con- 
stants in this model. From equations ([1]) and ©, we can 
obtain H c2 (T, 6). 

The diffusivity D\ along the c axis is smaller than the 
in-plane diffusivity D1 b in MgB2 due to the nearly 2D 
nature of the a band. On the other hand, the values of 
£>2 an( i ^ 2 b do not differ substantially because of the 
isotropic 3D nature of the ir band. The resulting rela- 
tions among diffusivities are D\ < Df and D% Df, 
which leads to the anomalous behavior of the anisotropy 
of H c2 (T). The in-plane diffusivity ratio Df/Df is an 
important parameter in the equation (TT|). 
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FIG. 3: (Color online) (a) Numerical ht to H C 2(T) obtained 
from transport data. The inset shows the temperature depen- 
dence of the anisotropy of H C 2- (b) Numerical fit to H C 2(9) at 
22 K and 32 K with the same parameters used to fit H C 2(T). 
From the fit, the diffusivity values of Di b =2.36 cm 2 /s and 
D% b =19.7 cm 2 /s were obtained; the coupling parameters ob- 
tained from the fit are L aa ~ 0.810, L ww ~ 0.285, L^^ ~ 0.25, 
and Laa ~ 0.18, respectively, close to the values obtained 
from ab-initio calculations (Ref. 34). The uncertainty of the 
fit parameters is no more than 5%, which is smaller than our 
experimental errors of 10%. 



We performed a numerical fit to three sets of trans- 
port data such as H c2 {T) at 9 =0°, H c2 (6) at T = 22 K, 
and T = 32 K using the equations (1) and (2), shown in 
Fig. [31 The relation between the best fit intraband diffu- 
sivities in the a and ir bands is = 8.5 x D\. This large 
i] = 8.5 is consistent with the absence of a sharp upward 



curvature in ifjj (T) at low T (see Fig. 1 in Ref. [2]), fre- 
quently observed in C-dopcd MgB2 with extremely high 
H c2 . The inset of Fig. [3ta) shows the anisotropy */h c2 (T) 
as a function of T. Again, this behavior is qualitatively 
consistent with that expected for 77 S> 1, see Fig. 3(c) in 
Ref. [2]. The superconducting critical field, H'^iO), for 
field applied parallel to the c axis, obtained from the fit, 
equals 10 T. This indicates the presence of strong multi- 
ple intraband scattering channels. The value of in-plane 
intraband diffusivity ratio r\ = 8.5 provides information 
about the type of the intraband scattcrers. The larger 
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FIG. 4: (Color online) The calculated thermodynamic critical 
field H c from the gap equations for two-band superconduc- 
tors. The inset shows the calculated gap values from the two 
band model. 



value of 77, (smaller value of Df) indicates the weakening 
of the 2D a band by certain types of impurities, such as 
C and N. These impurities affect the 2D landscape by 
replacing p xy orbitals of boron, and making the system 
more isotropic. The large value of DEJ compared to D" is 
also in good agreement with results from the a model^ 
and is the result of a large contribution of the n band to 
the total density of states. 



C. A and H c from the two-band model 

We calculated A using the parameters obtained from 
the H c2 (T,8) fit and the band calculations. The Lon- 
don equation for a two-gap superconductor is given by 



V x (A 2 V x H) 



tration depth is A£ 



H 



0, where the London pene- 



2 (T) = ne 2 fi (N 1 Df Aitanh^ + 
N^D^^taiah^): Indices of 1 and 2 represent the a 
band and the ir band, respectively. Ni and N2 are the 
electron densities of states. Ai and A2 are the gap mag- 
nitudes. Df b and D% b are the intraband diffusiviticsj 2 
Using Df b =2.36 cm 2 /s, Df = 19.7 cm 2 /s, obtained 
from the fit of H c2 (T,9), Ai(0) = 84 K, A 2 (0) = 33 K, 
obtained from the gap Eqs. (3), Ni = 0.3 states/a 3 eV, 
and N 2 = 0.41 states/a 3 eV with a unit cell volume of a 3 = 
87.2 A 3 — (obtained from the band calculations)^ we ob- 
tain Al(0) =170±10 nm, consistent with the measured 
value of A ab (0) =200±30 nm. The calculated X L (T) is 
shown as the blue curve in Fig. EJb) , consistent with the 
MFM experiment. This indicates that the two indepen- 
dent measurements of A(T) (MFM) and H c2 (T) (trans- 
port) in MgB 2 are complementary for investigating su- 
perconducting properties. 

The thermodynamic critical field (H c ) in MgB 2 is 



important for technological applications^ We evaluate 
H c using the band coupling parameters, obtained from 
H c2 (T,9), and the electron density of states obtained 
from the band calculations. 

The gap equations for two-band superconductors 3 ^- 



Ai 
A 2 



A^ 1 (0)A 1 y(A 1 ) 
7V 2 (0)A 2 F(A 2 ) 



0. 



■ are 



(3) 



: tanh 



2k B T 



with Y(Aj) = C" d£ . -- 1 

where g is the superconducting coupling matrix with 
.911 = NxL 22 /w,gx 2 = g 2 \ = NiL 12 /w = N 2 L 21 jw, 
and g 22 = N 2 Ln/w. lo c is some unknown cutoff fre- 
quency obtained from Eqs. (3) using the T c obtained 
from the transport data. Using the parameters obtained 
from the fit of H c2 (T, 9) , we have the superconducting 

0.46 -0.40 
-0.40 1.78 



coupling matrix, g = 



/(a 3 eV). The 



as 



free energy is calculated^ 

'cosh^i^^/TAl 



T = Y- 



A,5ijA*) - 
Then H r is 



given by H^/8n = -T . We calculate A X (T) and A 2 (T) 
as shown in Fig. 0] (inset). The calculated gap values 
at zero temperature are Ai(0) = 84 K and Ai(0) = 33 
K, which are slightly larger than reported values.— The 
thermodynamic critical field at zero temperature, calcu- 
lated from the two-band model, is approximately 2000 
Oe. This value is smaller than those previously obtained 
in polycrystalline MgB 2 by specific heat measurements 3 ^ 
and the values reported in clean single crystalsi 37 i 38 

As discussed earlier, the superconducting properties in 
multiband superconductors are affected by the interac- 
tions among the bands^i We obtain £ab(0) = 5.7 nm 
using H c2 (0) = $o/27r£ 2 (0) and our experimental value 

ffj 2 (0) = 10 T. We can then use the Ginzburg-Landau 
theory to estimate the thermodynamic critical field in 
the film, H c = $ /2V2ttA(0)£(0) = 2100 ± 300 Oe. This 
value is close to the calculated value of H c = 2000 Oe 
from the two band model. This suggests that the strong 
intraband scattering in the 3D 7r band makes the system 
more isotropic, and thus the system shows single band 
characteristics. 



IV. CONCLUSION 

In conclusion, we have measured X a b(T) and H c2 (T, 9) 
in a MgB 2 film. Our analysis of H c2 (T, 9) shows that the 
large value of the in-plane intra-band diffusivity in the 
3D 7r band is due to the presence of non-magnetic im- 
purities such as C and N, indicating the system is more 
isotropic, which is partly responsible for a large A. We 
calculated A and H c employing the gap equations for the 
two-band superconductors using the parameters obtained 
from H c2 (T, 9) and derived from band calculations. The 
calculated Al(0) =170±10 nm is close to the measured 
A(0) = 200 ± 30 nm in MgB 2 film, indicating that two 
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independent measurements, such as MFM and transport, 
are complementary, and provides a venue for thoroughly 
investigating superconducting properties. The determi- 
nation of H C (T) in clean MgB2 and in multi-band super- 
conductors in general is a fascinating problem with both 
fundamental and technological relevance. 
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